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Abstract 

The statistical theory of certain complex wave interference phenomena, like the 
statistical fluctuations of transmission and reflection of waves, is of considerable 
interest in many fields of physics. In this article we shall be mainly interested in 
those situations where the complexity derives from the quenched randomness of 
scattering potentials, as in the case of disordered conductors, or, more in general, 
disordered waveguides. 

In studies performed in such systems one has found remarkable statistical regulari- 
ties, in the sense that the probability distribution for various macroscopic quantities 
involves a rather small number of relevant physical parameters, while the rest of the 
microscopic details serves as mere "scaffolding" . We shall review past work in which 
this feature was captured following a maximum-entropy approach, as well as later 
studies in which the existence of a limiting distribution, in the sense of a gener- 
alized central-limit theorem, has been actually demonstrated. We then describe a 
microscopic potential model that was developed recently, which gives rise to a fur- 
ther generalization of the central-limit theorem and thus to a limiting macroscopic 
statistics. 
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1 INTRODUCTION 



Complex scattering of waves has captured the interest of physicists for a long 
time[l]. For instance, the problem of coherent multiple scattering of waves, 
which has long been of great importance in optics, has seen a revived interest 
in relation to the phenomenon of localization. 

The present article fits in the general topic of "Statistical theory of complex 
wave-interference phenomena". In particular, we shall study the statistical 
fluctuations of transmission and reflection of waves, which are of considerable 
interest in mesoscopic physics. Complexity in wave scattering may derive from: 

(i) the chaotic nature of the underlying classical dynamics, as in microwave 
cavities and quantum dots, or 

(ii) the randomness of the scattering potentials in a disordered medium, as a 
disordered conductor, or a disordered waveguide carrying classical waves 
(electromagnetic, elastic, etc.), on which we shall concentrate here. 

Why do statistics on the results of a scattering process? The point is that the 
interference pattern resulting from the coherent multiple scattering of waves 
from the systems described above is so complex (a small variation in some 
external parameter changes it completely) that only a statistical treatment is 
meaningful. 

We shall find a recurrent theme in our presentation: the statistical regularity 
of the behavior, which involves a relatively small number of relevant physical 
parameters, while the rest of details serves as mere "scaffolding". This feature 
was captured in the past following a Maximum-Entropy Approach, within 
the powerful, non-perturbative, framework known as Random-Matrix Theory: 
Shannon's information entropy is maximized, subject to the symmetries and 
constraints that are physically relevant [2,3]. Later, generalized Central-Limit 
Theorems (CLT) have been demonstrated [4,5]. Here we shall revisit past and 
recent efforts towards discovering universal features in the statistical scattering 
of waves in disordered waveguides. 

The paper is organized as follows. We first indicate the various physical regimes 
to be encountered in the problem of disordered conductors. We then mention 
how the statistical regularities in the problem have been captured in the past 
within a Maximum- Entropy Approach: we briefiy mention the transfer-matrix 
method that was used and the Random-Matrix Theory model that was con- 
structed, giving rise to a diffusion equation in transfer-matrix space. We then 
indicate a CLT that was proved, thus showing that the Maximum-Entropy 
Approach captures the universal features found in the CLT. We then describe 
a microscopic potential model that was developed recently, which gives rise 
to a further generalization of the CLT and thus to a limiting macroscopic 
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statistics. We then give our conclusions. 



2 THE MAXIMUM-ENTROPY APPROACH 



The fundamental physical process occurring in the system under study, the 
disordered waveguide shown schematically in Fig. 1(a), is coherent multiple 
scattering of waves. The length of the disordered section of the waveguide is 
L, W is its width (the system being assumed two-dimensional) and N is the 
total number of running modes supported by the waveguide; the mean-free- 
path (mfp) is denoted by £. The system might as well be a disordered conductor 
with the same dimensions. In the transport of waves through the system one 
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Fig. 1. A disordered waveguide and the various physical regimes. 

encounters various physical regimes, indicated in Fig. 1(b): i) the ballistic 
regime, for < L < ii) the diffusive, or metallic regime, for £ < L < Nl] 
iii) the insulating regime, for Nl < L. 

The aim is to calculate the transport properties of waves through this system. 
One important quantity is, of course the conductance of the disordered quasi- 
one-dimensional (q-ld) system, which is given by Landauer's formula: 



N 
a,b=l 



(1) 

(2) 



where a, b denote the transverse (running) modes: a,b — 1, - ■ ■ ,N. Lan- 
dauer's relation allows calculating the conductance -a transport quantity- 
from the scattering properties of the system. In microwave systems one can 
actually measure the individual transmission coefficients T^b = \t, 
as T = Y.ab=i l^aftP, experimentally. 



as well 
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The main approaches to the problem have been [1,6]: i) Perturbation theory 
in the disordered potential; ii) Supersymmetry methods, giving rise to a non- 
linear sigma model and iii) Random-Matrix Theory (RMT) models of the 
scattering matrix S, or the transfer matrix M of the system. Here we shall 
concentrate on these latter non-perturbative models. 

One has observed remarkable statistical regularities, in the sense that the 
probability distribution for various macroscopic quantities involves a relatively 
small number of relevant physical parameters (essentially the mean free path 
i). Within a Random-Matrix Theory scheme, this feature was captured in the 
past following a Maximum-Entropy Approach [2] , which we now describe. 

To a waveguide of length L we assign the transfer matrix M" (see Fig. 2), 
which has the property that acting on the wave amplitudes a^^^ on the left 
gives the amplitudes a^'^^ on the right: 

M"aW=a(2). (3) 



For every configuration of disorder we have one transfer matrix M". If we 
assign a probability density pl{M") to our transfer matrices, what results is 
a Random-Matrix Theory of transfer matrices. 

To the waveguide of length L we now add a "Building Block" (BB) of thickness 
5L which is much shorter than L, but still contains many weak scatterers (the 
so called dense- weak-scattering limit (DWSL)), as shown schematically in Fig. 
2. 




The transfer matrix M' associated with the BB will be considered to be sta- 
tistically independent from M"; we write its probability density as psl{M')- 
The transfer matrices M" and M' are combined as 

M = M'M" (4) 



to obtain the total transfer matrix of the combined system, while the proba- 
bility density of the latter, pl+sl{M), is obtained from the individual ones by 
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the "convolution" 

Pl+sl{M) = J pl{{M')-'M) psl{M') dji{M') = pl^Psl , 



(5) 



where dn{M') is the invariant measure associated with the group of transfer 
matrices. Eq. (5) has the structure of the Smoluchowski equation in Brownian 
motion theory [7]. 

Expecting the results to be largely independent of the details of the BB, the 
distribution psl{M') for the BB was modelled, in Ref. [2], by maximizing the 
Shannon entropy 

S\p\^- j p{M') \np{M')dM' , (6) 

subject to the constraints 

j psL{M')dfi{M') ~ 1 (normalization condition) (7) 
1 ^ AT" 

a, 6=1 

£ being the elastic mean- free-path (mfp). The only physical information con- 
veyed hy psl{M') is the mfp £ . 

Following the procedure developed in the theory of Brownian motion, we can 
convert the integral equation (5) into a differential equation, with the result 
[2,3]: 
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A„(1 + A„)J(^)(A) 



dXa J(^)(A) 



(9) 



This is a diffusion equation in transfer-matrix space, known as the DMPK 
equation (after Dorokhov [8] and Mello, Pereyra and Kumar [2]), which gov- 
erns the "evolution" with the length L of the waveguide of the probability 
distribution w^^^X) of our transfer matrices. This equation has been written 
in terms of new variables as follows: i) s = L/i denotes the length in units 
of the mfp; ii) A^ (a = 1, ■ ■ ■ , N) arc "radial" variables in terms of which the 
conductance can be written as g = Y,a=i 1/(1 + -^0)) iii) have used the Jaco- 
bian Jg(A) = Yla<b l-^a — A^l^. Finally, Eq. (9) has to be solved with the initial 
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condition Wq^\X) — 5{X). The quantity (3 (— 1,2,4) denotes the universahty 
class of random- matrix theory [9] . 

The distribution P{g) of the conductance is known to evolve from a Gaussian 
(deep in the diffusive, metallic regime) to a log-normal distribution (deep in 
the localized, insulating regime) [6]. Although P{g) cannot be easily obtained 
algebraically from the above expressions, various approximations show that 
this behavior is well described by the DMPK equation (9). In the crossover 
regime, Ref. [10] found the main statistical properties of P{g) arising from the 
DMPK equation using a Monte Carlo calculation. The results are shown in 
Fig. 3, where P{g) is plotted (solid lines) for both [3 = 1 and 2 for different 
values of {g). The results of extensive q-ld tight-binding-model calculations 
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Fig. 3. Conductance distribution from the DMPK equation and from numerical 
simulations for bulk and surface disorder, as explained in the text (after Ref. [10]). 

(symbols), carried out for bulk disorder, are also shown. We see that the 
agreement between the two types of results is excellent. In contrast, the results 
for surface disorder and (3 = 1, also shown in the figure for comparison, are 
not described properly by the DMPK equation. 



3 CENTRAL-LIMIT THEOREMS 



It was shown later [4] that a limiting distribution for ^^^(A) arises when the 
individual, microscopic, scattering units arc combined in the so called dense- 
weak-scattering limit (DWSL) and within a particular class of models. The 
DWSL corresponds to a large density of weak scatterers, with a fixed mfp i. 
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In the particular model that was studied, the individual scattering units were 
defined through their transfer matrices Mj and an "isotropic" distribution of 
their phases was assumed. The limiting distribution that arises depends only 
on £ and is insensitive to other details of the microscopic distribution: it thus 
constitutes a generalized central-limit theorem (CLT). The result turns out to 
be identical to the DMPK equation (9) found in the maximum-entropy model 
described above. We can thus say that the maximum- entropy model selects the 
limiting distribution, in the sense of the DWSL, within a class of models for 
the transfer matrices of the individual scattering units. 

A class of limiting distributions wider than that of Ref. [4] was studied by one 
of the present authors (PAM) and S. Tomsovic in Ref. [5], in which the isotropy 
assumption of Ref. [4] was relaxed to a large extent. In Ref. [5] the DWSL plays 
again an essential role and the result is a more general CLT than that of Ref. 
[4]. The evolution with L is described by a generalized diffusion equation, in 
which the diffusion coefficients are the inverse mfp's for the various scattering 
processes that may occur in the problem. When the various mfp's can be 
represented by a single one, one encounters the DMPK equation that was 
described above. Thus the model of Ref. [5] appears as a possible candidate to 
study the influence of the specific scattering properties of the various modes, 
which seem to be relevant, for instance, for the problem of waveguides with 
surface disorder, where DMPK does not give a proper description. 



4 A POTENTIAL MODEL FOR THE EVOLUTION OF EXPEC- 
TATION VALUES. A CENTRAL-LIMIT THEOREM 



We first present a general way of expressing the expectation value of an ob- 
servable when we add a BB to an already existing waveguide of length L, as 
shown in Fig. 2. The transfer matrix of the two pieces is combined as in Eq. 
(4), which can also be written as 

M = M" + = M" + sM" , (10) 

where we have expressed the transfer matrix of the BB as 

M' = 7 + £. (11) 



Consider now a function F{M) of the transfer matrix M, whose statistical 

properties we want to study: it might be, for instance, the conductance G 
studied earlier, the transmission coefficient Tab, or any other quantity of phys- 
ical interest. Its average for the enlarged piece can be written in terms of that 
for the original one as 
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{m))L+SL 



1^ d^F{M)\ 
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To simplify the notation, we have indicated symbohcally with dots the relevant 
summations over channel and block indices (see Ref. [5] for more details). 

It is clear that we now need an expression for the various moments of the 
quantities e associated with the BB. Recently [11], such an expression has 
been obtained from a potential model, which we now outline. We construct 
the BB as a sequence of m » 1 random 5-potential shoes, such that 



d -C 5L -C {A, 



(13) 



as shown schematically in Fig. 4. 
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Fig. 4. Construction of the BB using (5-potential slices in the regime defined by the 
inequaUties (13). 



The r-th 5-slice potential, defined as 
Ur{x,y) = Ur{y)S{x - Xr), 



(14) 



has matrix elements with respect to channels given by 



(15) 



in terms of which we specify the statistical model. The m potentials Ur, r ~ 
l,---,m, are assumed to be statistically independent, identically distributed, 
with zero average and, for simplicity, zero odd moments, so that 

{{UrU) = (16) 
((^irl)al6l(^ir2)a262) = «2(ai&l, a2&2)(^rir2 (17) 
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(('^n)ai6i('^r-2)a262('^r-3)a363('^r4)a464) — ^^2(0161, 02^*2) '«2( 03^3, Ci4b4)5rir2^r3r4 
+K2{aibi, azhz)K2{a2h2, aih^jSr^rz^riri 

+K2{a]hi, 0464) K2 (02^2, a'zbz)5r^rA^r2r3 

(18) 

' ' ' ) 

where ^2(0161, 0262), 't4(ai&i, 02^2, os&s, 04^4), etc., denote the second, fourth, 
etc., cumulants of {ur)ah- Prom these expressions we can calculate, in the 
DWSL, the various moments of e needed in Eq. (12). One finds that the first 
moment vanishes, the second moment behaves linearly with 5L and higher 
moments behave as higher powers thereof. Also, the very important result 
emerges that the dependence on the cumulants of the potential higher than the 
second drops out in the DWSL. The diffusion coefficients D^Jlj^, or inverse 
mfp's lahi defined as 

("«^"cd)^,5L = 2<c^(^'^)^^ + ---' (19) 



depend only upon the second cumulants of the potential. They are energy 
dependent and also length dependent. 

Finally, we take the first term on the r.h.s. of Eq. (12) to the l.h.s., divide 
both sides hy 5L and take the limit (5L — > 0. The result is the Fokker-Planck 
equation: 

d{F{M))^ _ jk,lm/ i\ jkn n jmp d'^F \ 



The fact that cumulants of the potential higher than the second are irrelevant 
in the end signals the existence of a generalized CLT: once the mfp's are 
specified, the limiting equation (20) is universal, i.e., independent of other 
details of the microscopic statistics. 

One of the main difficulties in solving Eq. (20), both analytically and nu- 
merically, is that it involves averages of different quantities on the l.h.s. and 
on the r.h.s. So far, that equation has been solved analytically for the one- 
open-channel case {N = 1) and a restricted number of "observables" only: 
the results thus obtained are in excellent agreement with microscopic calcu- 
lations [11]. Numerically, we have found no "direct" way of solving Eq. (20). 
Recently, a numerical algorithm which was called "random walk in transfer- 
matrix space" has been implemented [11,12], in which a BB is constructed 
with the property given in Eq. (19), and then combined with successive BB's 
to construct a waveguide of finite length L. The results have been compared 
with those arising from microscopic calculations, in which the entities that are 
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combined are, literally, individual potential slices. For situations in which we 
have bulk disorder, the comparison is excellent, even for quantities which are 
not described properly by DMPK. For surface disorder, preliminary results 
indicate a reasonable agreement for the observables examined so far. 



5 CONCLUSIONS 

In this paper we have first revisited some earlier results which seem to indicate 
that a maximum-entropy approach to the problem of transport in disordered 
waveguides works well when there is a central-limit theorem (CLT) "behind 
the scenes" . 

We have briefly described a rather general CLT that was obtained in the past, 
in Ref. [5]. 

Recently, we have shown that a CLT arises in a model consisting of a ran- 
dom distribution of (5-potential slices. The parameters on which the result 
depends are the mfp's -which depend on the variances of the potential- 
slices matrix elements- and the correlation coefficient of these matrix elements. 
Other details of the potential distribution arc "washed out" in the dense-weak- 
scattering limit (DWSL). The result is expressed in terms of a generalized dif- 
fusion equation for the evolution with length of expectation values of physical 
observables. 

Numerical results based on the "random walk in transfer-matrix space" method 
indicate an excellent agreement with microscopic calculations with bulk dis- 
order, and a reasonable one for the problem with surface disorder. 

More effort is needed towards an analytical, as well as a numerical, treatment 
of the diffusion equation. 
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